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Abstract 

Denote by w(A) the numerical radius of a bounded linear operator A acting on Hilbert space. 
Suppose that A is invertible and that w(A) < 1+e and w(A^ 1 ) < 1+e for some e > 0. It is shown 



(N 



that inf{||A-Z7|| : U unitary} < ce 1 / 4 for some constant c > 0. This generalizes a result due to 
J.G. Stampfli, which is obtained for e = 0. An example is given showing that the exponent 1/4 is 
optimal. The more general case of the operator p-radius w p {-) is discussed for 1 < p < 2. 

1 Introduction and statement of the results 

< 

Ph " Let H be a complex Hilbert space endowed with the inner product (•, •) and the associated norm || • ||. 

We denote by 13(H) the C*-algebra of all bounded linear operators on H equipped with the operator 

> ■ 

, norm 

\\A\\ = su V {\\Ah\\ : h G H,\\h\\ = 1}. 
It is easy to see that unitary operators can be characterized as invertible contractions with contractive 



> 

m 
o 



inverses, i.e. as operators A with \\A\\ < 1 and \\A 1 || < 1. More generally, if A G B(H) is invertible 
then 



^ . inf : U unitary } = max \\A\\ - 1, 1 - t— - 



A- l \ 

We refer to [6] Theorem 1.3] and [9l Theorem 1] for a proof of this equality using the polar decompo- 
sition of bounded operators. It also follows from this proof that if A G B(H) is an invertible operator 
satisfying \\A\\ < r and ||^4 _1 || < r for some r > 1, then there exists a unitary operator U G B(H) 
such that 1 1 -A— < r—1. 

The numerical radius of the operator A is defined by 

w(A) =sup{\(Ah,h)\ : h G H, \\h\\ = 1}. 

Stampfli has proved in [8] that numerical radius contractivity of A and of its inverse A -1 , that is 
w(A) < 1 and w(A^ 1 ) < 1, imply that A is unitary. We define a function ip(r) for r > 1 by 

^(r) = sup{||A|| : A G B{H),w{A) < r, wiA" 1 ) < r}, 

the supremum being also considered over all Hilbert spaces H. Then the conditions w(A) < r and 
w(A^ 1 ) < r imply max(||A|| — 1, 1— ||A _1 || _1 ) < max(||A||— 1, ||A _1 || — 1) < ip(r)— 1, hence the existence 
of a unitary operator U such that \\A— U\\ < ip(r)—l. We have the two-sided estimate 



r+V r 2 — 1 < ^/>(r) < 2r. 

The upper bound follows from the well-known inequalities w(A) < \\A\\ < 2w(A), while the lower 
bound is obtained by choosing H = C 2 and 



A= (p -l) With y = v/r2 ~ 1 ' 
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in the definition of ip. Indeed, we have A = A 1 , w(A) = \fl + y 2 = r, and ||A|| = y + \/l+y 2 
r + \/r 2 — 1. 

Our first aim is to improve the upper estimate. 
Theorem 1.1. Let r > 1. Then 



il){r) < X(r) + y/X{r) 2 -l, with X(r) = r + \J r 2 - 1. (1) 

The estimate given in Theorem II .11 is more accurate than ifj(r) < 2r for r close to 1, more precisely 
for 1 < r < 1.0290855... . It also gives = 1 (leading to Stampfli's result) and the following 

asymptotic estimate. 

Corollary 1.2. We have 

ip(l+e) < 1 + ^8e + 0{e 1/2 ), e -> 0. 

Our second aim is to prove that the exponent 1/4 in Corollary [L2] is optimal. This is a consequence 
of the following result. 

Theorem 1.3. Let n be a positive integer of the form n = 8k + 4. There exists a n x n invertible 
matrix A n with complex entries such that 

w{A n )<-^—, w ( A - l )<-^, \\A n \\ = \+ ' 



cos £ cos 2 8Vn 

Indeed, Theorem 11.31 implies that 

^)>Pn|| = l + - 



Taking 1+e = 1/ cos \ = 1 + ^ + O(^), we see that the exponent \ cannot be improved. 

More generally, we can consider for p > 1 the p-radius w p (A) introduced by Sz.-Nagy and Foia§ 
(see (5j Chapter 1] and the references therein). Consider the class C p of operators T G B(H) which 
admit unitary /^-dilations, i.e. there exist a super-space % D H and a unitary operator U G B{%) 
such that 

T n = pPU n P*, for n = 1,2,.... 

Here P denotes the orthogonal projection from Ti onto H . Then the operator p-radius is defined by 

w p (A) = inf{A > ; \~ 1 A G C p }. 

From this definition it is easily seen that r(A) < Wp(-A) < where r(A) denotes the spectral 

radius of A. Also, w p (A) is a non-increasing function of p. Another equivalent definition follows from 
[51 Theorem 11.1]: 



w p (A) = sup {(1-1) \(Ah,h}\ + ^i-l)2\(Ah,h)\ 2 + ( 2 p -l)\\Ah\\ 2 } : 
E p = {h G H ; = 1 and(l-I) 2 |(A/ l ,/ l )| 2 - (1-|) ||^|| 2 > 0}. 



with 



P I ) II "II V p/ I \ l/l V p ) 

Notice that £ p = {h G H ; \\h\\ = 1} whenever 1 < p < 2. This shows that wi(A) = \\A\\, W2(A) = w(A) 
and Wp(A) is a convex function of ^4 if 1 < p < 2. 
We now define a function ip p (r) for r > 1 by 

^p(r) = sup{\\A\\ ; A G B(H),w p (A) < r, ^(A" 1 ) < r}. 

As before, the conditions u>p(A) < r and ^(A -1 ) < r imply the existence of a unitary operator {7 
such that \\A— U\\ < ip p {r)— 1, and we have tp p (r) < pr. We will generalize the estimate ([T|) from 
Theorem 11.11 by proving, for 1 < p < 2, the following result. 
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Theorem 1.4. For 1 < p < 2 we have 

Mr) < X p {r) + ^X p {rf - 1, (2) 

2 + pr 2 - p + y^ + pr 2 -p) 2 -4r 2 

with Xp{r) = v -— . 

2r 

Corollary 1.5. For \ < p <2 we have 

^>(l+e) < 1 + \/8(p - l)e + 0{s l l 2 ), e ->■ 0. 

We recover in this way for 1 < p < 2 the recent result of Ando and Li [21 Theorem 2.3], namely 
that w p (A) < 1 and w p (A~ l ) < 1 imply A is unitary. The range 1 < p < 2 coincides with the range 
of those p > 1 for which w p {-) is a norm. Contrarily to [2\, we have not been able to treat the case 
p > 2. 

The organization of the paper is as follows. In Section 2 we prove Theorem 11.41 which reduces to 
Theorem II .11 in the case p = 2. The proof of Theorem 11.31 which shows the optimality of the exponent 
1/4 in Corollary 11.21 is given in Section 3. 

As a concluding remark, we would like to mention that the present developments have been in- 
fluenced by the recent work of Sano/Uchiyama [7j and Ando/Li [2J. In [3], inspired by the paper of 
Stampfli [8], we have developed another (more complicated) approach in the case p = 2. 



2 Proof of Theorem 11.41 about ift p 

Let us consider M = \{A + (A*)~ l ); then 

M*M - 1 = \{A*A + (A* A)' 1 - 2) > 0. 

This implies ||M _1 || < 1. In what follows C 1 / 2 will denote the positive square root of the self-adjoint 
positive operator C. The relation (A* A - 2M*M + l) 2 = 4M*M(M*M - 1) yields 

A* A - 2M*M + 1 < 2(M*M) 1 / 2 (M*M - I) 1 ' 2 , 

whence A* A < ((M*M) 1/2 + (M*M - 1) 1/2 ) 2 . 

Therefore \\A\\ < ||M|| + ^Mll 2 - 1. 

We now assume 1 < p < 2. Then w p (.) is a norm and the two conditions w p {A) < r and 
w p {A~ l ) < r imply w p (M) < r. The desired estimate of ip p (r) will follow from the following auxiliary 
result. 

Lemma 2.1. Assume p > 1. Then the assumptions w p (M) < r and ||M _1 || < 1 imply \\M\\ < X p (r). 

Proof. The contractivity of M _1 implies 

|M| < ||M«||, (yueH). (3) 

As w p (M) < r, it follows from a generalization by Durszt [1] of a decomposition due to Ando |1], that 
the operator M can be decomposed as 

M = prB 1/2 UC 1/2 , 

with U unitary, C selfadjoint satisfying < C < 1, and B = /(C) with f(x) = (1— x)/(l— p(2— p)x) _1 . 
Notice that / is a decreasing function on the segment [0,1] and an involution: /(/(x)) = x. Let 
[a, /3] be the smallest segment containing the spectrum of C. Then [^/a, \//3] is the smallest segment 
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containing the spectrum of C 1//2 and [y/JJfi), \J f(ct)] is the smallest segment containing the spectrum 
of B 1 / 2 . We have 



\\u\\ < \\Mu\\ < pr y/f(a)\\C 1/2 u\\, (V« 6 H). 

Choosing a sequence u n of norm-one vectors (||itn|| = 1) such that ||C 1//2 u n || tends to y/a, we first get 
1 < pry / af(a), i.e. 1 — (2+pr 2 —p)pa + p 2 r 2 a 2 < 0. Consequently we have 



2+pr 2 -p- J(2+pr 2 -p) 2 -4r 2 2+pr 2 -p+ J (2+pr 2 -p) 2 -4r 2 

, 5 — a — , o • 

2pr 2 2pr l 

and by a = /(/(a)) 



2+pr 2 -p - v / (2+pr 2 -p) 2 -4r 2 2+pr 2 -p + y/(2+pr 2 -p) 2 -4r 2 

2p~7 2 " /(a) " 2p~7 2 ' 

Similarly, noticing that || (Af*)- 1 1| < 1, M* = prC l ' 2 U*B 1 / 2 and C = f(B), we obtain 



2+pr 2 -p - yj(2+pr 2 -p) 2 -4r 2 < ^ < 2+pr 2 -p + y/ ' (2+pr 2 -p) 2 -4r 2 
2 pr 2 ~ — 2 pr 2 

Therefore 



||M|| < pr \\ B W\\ \\CV 2 \\ = pr < 2 +P±R±V^E^E^ , 

This shows that \\M\\ < X p (r). □ 



3 The exponent 1/4 is optimal (Proof of Theorem LL3J) 

Consider the family of n x n matrices A = DBD, defined for n = 8k + 4, by 

D = diag(e i7r / 2n , . . . , e ( 2 ^W2n . . . ; e (2n-l)W2n ); 

B = I + n \/o E, where E is a matrix whose entries are defined as 
eij = 1 if 3A; + 2 < \i — j\ < 5k + 2, = otherwise. 

We first remark that \\A\\ = \\B\\ = 1 + ^=. Indeed, B is a symmetric matrix with non negative 
entries, Be = (l+g~^) e with e T = (1, 1, 1 . . . , 1). Thus \\B\\ = r(B) = 1+-^= by the Perron-Probenius 
theorem. 

Consider now the permutation matrix P defined by pij = 1 if z = j + 1 modulo n and p$j = 
otherwise and the diagonal matrix A =diag(l, . . . , 1, —1). Then P~ 1 DP = e l7T ^ n AD and P~ 1 EP = 
E, whence (PA^APA = e 2i7T / n A. Since PA is a unitary matrix, the numerical range PU(A) = 
{{Au,u),;\\u\\ = 1} of A satisfies W(A) = ^((PA^APA) = e 2ilx l n W(A). This shows that the 
numerical range of A is invariant by the rotation of angle 2ir/n centered in 0, and the same property 
also holds for the numerical range of A . 

We postpone the proof of the estimates |||^4 + A*)|| < 1 and || \{A~ l + (yl _1 )*)|| < 1 to later 
sections. Using these estimates, we obtain that the numerical range is contained in the half- 

plane {z ; Re z < 1}, whence in the regular n-sided polygon given by the intersection of the half-planes 
{z ; Re(e 2mk l n z) < 1}, k = l,...,n. Consequently w(A) < l/cos(7r/n). The proof of w(A^ 1 ) < 
l/cos(7r/n) is similar. 
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3.1 Proof of < 1. 

Since the (^,j)-entry of A is e^ + i~^ l ~ (dtj + 2 ^3% ) ' ^ e ma ^ r i x \{A + A*) is a real symmetric 



matrix whose entry is cos ((i+j— 1)^) y5ij + 2 ^3/ 2 j • It suffices to show that, for every u = 

(ui, • • • ,u n ) T G R n , we have ||u|| 2 - Re(Au,u) > 0. Let £ = ;1 < i,j < n,3k + 2 < \i - j\ < 

5k + 2}. The inequality which has to be proved is equivalent to 

n 

^2sin 2 ((i-I)^K 2 -^ J cos((i+j-l)^)u^ > 0. 

Setting = Uj sin((j — ^)^), this may be also written as follows 

2|M| 2 - (Mv, v} + ^ (Ev, v) > 0, (v G K n ). (4) 
Here M is the matrix whose entries are defined by 



m 



j/2 cot ((*-5)n) cot ((i-^)^)' if (*>J) G m u = otherwise. 



% 3 ~ 2n 

We will see that the Frobenius (or Hilbert-Schmidt) norm of M satisfies \\M\\ F < ^9/32 < 3/4. 
A fortiori, the operator norm of M satisfies \\M\\ < |. Together with \\E\\ = n/4, this shows that 
11-^11 + 3/2 11-^11 — I- Property (j4|) is now verified. 

It remains to show that ||M|| 2 < First we notice that m,j = mji = m n+ i-i jn+ i-j, and ma = 0. 
Hence, with £' = G £ ;i < j and i + j < n + 1}, 

||M|| 2 =2^|m u -| 2 <4 ^ K/. 

We have, for (i, j) G 

2j < « + j + 5A; + 2 < n + 5/c + 3 = 13k + 7, thus 3/c + 3<j<±^±?, 
2«<i + j-3/c-2<n-3/c-l = 5A; + 3, thus 1 < i < ^±2. 

This shows that 

3?r ^ 3fc+2 < / • 1\7T <- 13fc+6 < 3?r V| PT1rP I CCl + ((n — I~l I < C nt ^ < ^ 

16 — 16^+8^ — V 2^n - lefc+S 7 ' - 77 16' rlence I co HU 2 J n ^ - LQl 16 - 2 ■ 

We also use the estimate cot((i— 3)-) < n/(ir(i— 5)) and the classical relation £)i>i(* — 1/2) -2 = 7r 2 /2 
to obtain 



IMII 2 < 4 V |m i? -| 2 < -A? ^ V 7 l — -k (2k+l) 9 

(ij)e£' i> 



4n 3 ^2 (i- 1/2) 2 v ^ y 4 32' 



3.2 Proof of + < 1. 

We start from 

(A- 1 )* = D(l + ^ m E)- 1 D 

= D 2 - ^DED + ^DE\l + ^E)-*D, 

and we want to show that ||it|| 2 — Re(A~ 1 u,u) > 0. As previously, we set vj = Uj sin((j — i)^). The 
inequality + < 1 is equivalent to 

2|H| 2 - ((Mi + M 2 + M 3 + Af 4 )«, v) > 0, (u G R n ). 
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Here the entries of the matrices M p , 1 < p < 4, are given by 
(mi)ij = -2^72 cot((i-i)f) cot((j-i)^) eii , 

{ m 2)ij = 2n 3/2 e «J) 

(^3)ij = 47? «>t((i-£)£) COt(0-^)S)Ai. 
(777.4) ij = —J^Z fij > 

eij and fij respectively denoting the entries of the matrices E and F = E 2 (l + j^ i?) -1 . Noticing 

that Mi = — M, we have HMJ < f , ||M 2 || = ^, ||F|| < < fj and ||M 4 || = ^ ||F||. Now 

we use 

||M 3 || 2 < ||M 3 || 2 F < ^nu«|/ f X;|oot((t-|)S)| 2 |oot(0--^)| 2 , 
together with 

E 1 cot ((*- ^i 2 1 oot(o - \)i )i 2 = ( E 1 ^-\)i)\ 2 ) 2 

i,j i=l 

«/2 2 
<4(E|cot((7-i)^)| 2 ) < 77 4 , 
1=1 

to obtain 

||M 3 ||<^max|^|. 

Using the notation H^Hoo := max{||.E'ii|| 00 ; w G C n , ||u||oo < 1} for the operator norm induced by the 
maximum norm in C rf , it holds H-E^oo = 77/4, whence || 2ra 3/2 ^lloo < 1/8 and thus ||(1 + ^E)' 1 ^ < 
1 _ 1 1 y 8 = |. This shows that 

In 

maxl/^l < ||(1 + ^TjSj-^looinaxle?-! < y, 

by denoting e 2 - the entries of the matrix E 2 and noticing that maxjj |e? | = n/4. Finally, we obtain 
||M 3 || < ± and HM1 + M2 + M3 + M4H + ! + £ + ^ < 1. 
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